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Abstract 



Attention has been brought to the possibility that statistical fluctuation 



' properties of several complex spectra, or, well-known number sequences may 

■ display strong signatures that the Hamiltonian yielding them as eigenvalues is 

PT-symmetric (Pseudo-Hermitian). We find that the random matrix theory 



of pseudo-Hermitian Hamiltonians gives rise to new universalities of level- 
spacing distributions other than those of GOE, GUE and GSE of Wigner and 
Dyson. We call the new proposals as Gaussian Pseudo-Orthogonal Ensemble 
and Gaussian Pseudo-Unitary Ensemble. We are also led to speculate that 



the enigmatic Riemann-zeros (^ ± it n ) would rather correspond to some PT- 
symmetric (pseudo-Hermitian) Hamiltonian. 

I. INTRODUCTION 

A large body of spectra of the bound levels and resonances is available in nuclear physics 
wherein the nuclear interaction Hamiltonian is unknown. The well-known prime numbers 
2,3,5,7,9,11,13,17,19,... do not have a representation so far. One may wonder if there is a 
Hermitian Hamiltonian that can yield them as its discrete eigenvalues. The zeta function, 
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C(z), which is real on real line as per the hypothesis of Riemann (1859) has all its non-real 
zeros as \ ± it n {t n being real) [1]. The marathon and accurate computation of more than 
10 20 zeros of ((z) by Odlyzko [2] testifies RH the best, with not even a single exception so 
far. First few zeta-zeros (Riemann-zeros : RZs) are given by t\ = 14.13, ti = 21.02, t% = 
30.42, t 4 = 37.58. Hilbert and Polya have conjectured that t n could be like the eigenvalues 
of a Hermitian Hamiltonian. Consequently the completeness of the spectra will lead to the 
proof of one of the most enigmatic and formidable problems called RH. 

In nuclear spectroscopy, the pressing need was to identify the universality that under- 
lies the nuclear levels. The Nearest Neighbour Level Spacing Distribution (NNLDS) is the 
statistical distribution of the fluctuations around the mean level spacing of a collection of 
spectra under a class of fixed parity or other quantum numbers. Random Matrix Theory 
(RMT) was discovered in the late 1950s to predict universalities of NNLSD in various situa- 
tions. Among the notable names we have Wigner, Landau, Dyson, Gaudin, Mehta, Porter, 
Ginibre and Pandey to associate with RMT [3]. 

Since the interaction Hamiltonian is not known it would rather be taken as non- 
integrable, then there are three types of NNLSDs obtained by Wigner and Dyson. These are 
given as in Eq. (4) (see below) and called as GOE,GUE and GSE spacing-statistics [3]. Most 
remarkably the nuclear levels with same (J, 7r) are well-known to follow GOE statistics. The 
energy levels of the chaotic Sinai-billiard are known to follow the same statistics. This has 
supported the idea that chaos may be there in nuclear dynamics too. Sinai-billiard refers to 
a particle in a square region with hard, reflecting edges along with a hard, reflecting circle 
in its center. 

Well before the advent of RMT, Hilbert had prophesied that the real part of RZ are 
distributed as the eigenvalues of certain random Hermitian matrices. In fact, in terms of 
RMT it means that the real parts of RZs obey NNLSD corresponding to GUE. Over a 
hundred years old, this prophecy of Hilbert has been testified by Odlyzko [2] as late as 1989 
using more than 10 20 RZs. Montgomery in 1973 (see in Mehta [3]) analytically derived 
two-point correlation function of the real parts of RZs which turned out to be the same as 
that of GUE. Dyson had already expected this and it has also been confirmed numerically 
by Odlyzko [2]. These two remarkable affirmations have strengthened the Hilbert-Polya 
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conjecture to look for a Hermitian Hamiltonian for the RH. 

Last few years have witnessed an interesting phenomenon whereby the real discrete eigen- 
values need not necessarily be possessed by Hermitian Hamiltonians. Non-Hermitian, PT- 
symmetric [4] or pseudo- Hermitian [5,6] Hamiltonians too can possess real discrete spectrum. 
In RMT, the matrix ensembles GOE, GUE, GSE refer to Hamiltonians with TRI (Time Re- 
versal Invariance), without TRI, and with TRI including Kramer's degeneracy respectively. 
Recently, we have developed the Gaussian-random ensembles of pseudo-Hermitian matri- 
ces, that give rise to new "universalities" of NNLSD. We have called the new ensembles as 
GPUE [7,8] which are expected to represent the cases where Parity (P) and Time-reversal 
(T) symmetries are individually broken but preserved jointly. 

In this paper, we would like to present GPUE with more refinements and reorientations. 
In Section 2, we briefly introduce RMT with the ensembles of Wigner and Dyson. We then 
find a natural scope to go in for new ensembles. In Section 3 and 4 the new ensembles are 
described. In Section 5, we report two interesting dichotomies where spacing statistics are 
like GOE and GUE, despite Hamiltonians being pseudo-Hermitian. Following one of these 
dichotomies, in Section 6, we speculate on the possible features of Hamiltonian corresponding 
to RH. We present a summary of conclusions in Section 6. 

II. WIGNER-DYSON ENSEMBLES OF GAUSSIAN-RANDOM MATRICES : 

GOE, GUE, GSE 

An eigenspectrum is (practically) an outcome of the diagonalization of a Hamiltonian 
matrix. Since the analysis of NNLS requires at least two eigenlevels, in RMT, one begins 
with 2x2 Hamiltonian matrix for simplicity. The RMT takes an important note of the fact 
that for systems with TRI, the matrix Hamiltonians are real-symmetric (H R ); for systems 
without TRI, the Hamiltonians are Hermitian matrices (H H ) and for systems with TRI 
plus Kramer's degeneracy, the Hamiltonians are even (at least 4) dimensional (H K ). These 
Hamiltonian matrices are given as [3,9] : 
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a 7* -5 
a 5* 7 

(1) 

7 5 (3 
-5* 7* /3 j 

where a = a + b,j3 = a — b, r y = c + id,5 = e + if. In order to have a large collection of 
eigenvalue pairs, we assume the entries a, b, c, d, e, / are real and drawn independently from 
a Gaussian-random population. Notice that respective eigenvalues for (1) are 

E£ 2 = a±Vb 2 + c 2 , E* 2 = a±Vb 2 + c 2 + d 2 , E* 2 = a ± yjl? + c 2 + d 2 + e 2 + f 2 . (2) 

In RMT, to reemphasize, the energy eigenvalues for TRI systems turn out to have the form 
as E R , for non-TRI systems they are as E H and for TRI systems with Kramer's degeneracy 
are as E K s. The level-spacings (s = \Ei — E 2 \) in various cases are 

S R „ ^/WT^, s H ~ Vb 2 + c 2 + d 2 , s K ~ ^jb 2 + c 2 + d 2 + e 2 + f 2 . (3) 

Now the question to be asked is : what is the probability distribution of the level-spacing 
s, when a,b,...,f are Gaussian-random variables? More importantly one wants to know 
whether the levels have tendency to repel or attract each other and then what the degree of 
repulsion/attraction is. The real-symmetric matrices have an Orthogonal symmetry which 
corresponds to SO(N) : GOE. The Hermitian matrices have unitary symmetry which cor- 
responds to SU(N) : GUE. The matrices H K have Symplectic symmetry corresponding to 
Symplectic group S P (N) : GSE. These facts are used to derive Wigner-Dyson universalities 
in NNLSD as [3,9] 

2 QO 2 Ql8 

P°oe {x) = l X e-1-, P GUE (x) = ^x 2 e-^r, P GSE (x) = ^x 4 e"i^. (4) 

Here x is a scaled spacing (s) with respect to its mean value ((s)). For smaller values of 
spacing s, notice the tendency of level repulsion as P(x — > 0) — > 0. The degree of repulsion 
is linear, quadratic and quartic respectively for GOE, GUE and GSE. Most importantly 
these distributions (4) turn out to be excellent approximants to P(x) for N x iV matrices 
[3,9]. 

An alternative intuitive way of looking into these NNLSD is to see that the three statistics 
(4) correspond to the following multidimensional integral 
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R 



a + b c 
c a — b 



H 



H 



a 7 
7* P 



H 



K 



/oo p oo c oo 
/ / e - (a2+b2+c2 1 S(s - VWTd 2 ) da db dc, (5) 
-oo J — oo J — oo 

and similarly others for GOE and GSE. 

The common feature of the eigenvalues (2) or spacings (3) lies in their absolute reality. Let 
us ask the following questions : Firstly, can there be Hamiltonians possessing conditionally 
real eigenvalues or spacings e.g., s ~ \Jb 2 — c 2 (real iff b 2 > c 2 )and s ~ \fb 2 + c 2 — d 2 (real 
iff b 2 + c 2 > d 2 )? Secondly, what are the symmetries of such Hamiltonians? Thirdly, what 
are new universalities of NNLSD ? Such questions have led us to think of new Gaussian- 
random ensembles of pseudo- Hermit ian Hamiltonians [7,8]. By a Gaussian ensemble, we 
would mean that the probability distribution of Hamiltonian H is commonly given as 

P{H)=Ne b>— . (6) 



III. GAUSSIAN PSEUDO-ORTHOGONAL ENSEMBLES (GPOE) 



Pseudo-symmetric or complex-symmetric matrix Hamiltonians 



Let us consider a matrix Hamiltonian H given below which is pseudo-Hermit ian as r^Hr] 1 = 
if* [5] and pseudo-real i.e., pHp^ 1 = H*[ll]. It is self-pseudo-adjoint or symmetric as 
H' = H . Here r\ are p are preferably involutary operators. It has got conditionally real 
eigenvalues iff b 2 > c 2 Here prime, asterisk (/C ) and dagger denote transpose, conjugate and 
transpose-conjugate, respectively. 



H = 



a + b ic 
ic a — b 



, b 2 > c 2 , rj = p 



1 
-1 



a ± Vb 2 — c 2 . 



(7) 



One can construct an antilinear commutant 6 = p 1 /C [10] of H such that [0, H] = 
or QHQ- 1 = H and 2 = 1. We would like to assert that here P = p' 1 and T = JC 
and hence the antilinear symmetry 6 = PT. When eigenvalues are real (b 2 > c 2 ), we 
have PT^ n = (— l) n \l/ n . When b 2 < c 2 , the PT-symmetry is spontaneously broken. This 
Hamiltonian in our opinion is another realization of Hamiltonians with antilinear symmetry 
as visualized by Haake (page 217 in [9]) as [A,T>] = such that A 2 = 1. 
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We define pseudo-orthogonal transformation as p O' = 5 p O _1 5 _1 such that for any two 
arbitrary vectors from a linear space the scalar product remains invariant i.e., x'5y = x'Sy, 
where x = p Ox and y = p Oy. Let us represent p O, energy-eigenvalue matrix E and a 
metric 8 





cosh 9 i sinh 9 




-i 


P o = 




, s = 




—i sinh 9 cosh 9 




i 



E 



E l 
E 2 



(8) 



where — oo < 9 < oo. The single parameter matrix p O, is expressible as exp(2i9J 2 ) with 
■h = \iv y , constitutes a subgroup of S77(l,l) [12]. A very important consequence of the 
group connection is that we can generate all possible H in (7) as H = p \0 E p O _1 . This 
provides us with a unique connection between (a, b, c) and (Ei,E 2 ,9) and the consequent 
Jacobian is J = ^. We have 



a — — ^ — -, b= — — ; — - cosh 2^, c= ^— ; — ^■sinh26 ) . 



(9) 



2 2 2 

For brevity, we have used s = E x — E 2 and t = E x + E 2 . We can write the probability 
distribution, P{H) (6) for the Hamiltonian in (7) as 



P(a,b,c) = He- 



rn 



Using (9) we can transform (10) in terms of (t, s, 9), further integration over 9 on [— oo, oo] 
gives Joint Probability Distribution function of (Ei,E 2 ) as 

/ „2 \ ,2 



P{E 1 ,E 2 )=N'sK Q 



Aa 2 



e 4 CT ^ 



(11) 



Next the integration over t on [— oo, oo] yields the NNLSD as 

P(s) = N" s K 



2a 2 



(12) 



By finding (s) using (12) and introducing x = -^y, we eventually find the normalized NNLSD 
and call it as P GPOE {x) 



jGPOE, 



x) = 



327T 3 ° I 7T 2 



(13) 



If we write it as P GPOE (x) = a x K (px 2 ), we have a = 0.5818 and (5 = 0.4569. When 
< x < 0.5, we have P GPOE (x) ~ (0.5 — 1.2\nx)x. For any other pseudo-symmetric or 



complex symmetric matrix Hamiltonian that is composed of three independent Gaussian- 
random variables (a, b, c) appearing linearly in H, we claim that P GPOE (x) is the universality. 
The new "universality" shows a distinctly different behaviour as compared to the usual ones 
(see Fig. 1(b)). 



IV. GAUSSAIN PSEUDO-UNITARY ENSEMBLES 



Pseudo-Hermitian matrix Hamiltonians 



We now consider pseudo-Hermitian matrix Hamiltonians with four parameters (a, b, c, d) 



H = 



a + b d + ic 
-d + ic a — b 



^ e 2 = b 2 -c 2 + d 2 >0,r] = 



1 
-1 



, E h2 = a±e. 



(14) 



Here we have rjHi]' 1 = W , P and T operators can be constructed as prescribed in [10] an and 
antilinear commutant, 6, of H can be constructed as prescribed in [11]. Consider a trans- 
formation P U which preserves the pseudo-norm as x^rjy = x^rjy, where x — p Ux,y = p Uy. 
In doing so P U would satisfy an interesting condition i.e, P W = r\ p U~ 1 r]'~ 1 which is called 
pseudo-unit arity. (see e.g., [7]) 

A general three parameter (8, ip, <p) matrix, P U, which is pseudo-Unitary under the same 
metric i] (14) can be written as 



e 1 ^ cosh 9 e 1 ^ sinh 9 



,0 < < 2vr,0 < 9 < oo. 



(15) 



This constitutes a Lie group SU (1,1) [12] with generators as Jo = \&z, J\ = \i&y>Jz = 
— \i(j x . However, in order to construct the pseudo-Hermitian matrix (14) we require only 
two parameters in P U . The same situation arises [3,9] in case of GUE, where two out of three 
parameters suffice in writing the unitary matrix U , nevertheless it requires three-parameters 
to have SU (2). Thus, we take ip = in (15) and generate H in Eq. (14) as P UE p U~ l = H. 
This is how we go over to (E\, E 2 , (f>, ip) from (a, b, c, d). We find 

t s s s s 2 

a = -, 6=-cosh26 l , c— — sinh 29 cos 6, d = - sinh 29 sin 6 and J = — sinh 29. (16) 
2 2 2 2 4 
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The probability function (6) for H (14) works to P(a, b, c, d) = N exp[-(a 2 + b 2 +c 2 + d 2 )/a 2 ]. 
A similar procedure (yet more involved) as done in Section 3, from Eq. (10) to (13), leads 
us to a new NNLSD as 

P*™ (*) = * x erfc , B= g^i^. (17) 

If we write as P GPl/£ (:r) = a x e^ 2 erfc (7a;) where a = 2.5433, /3 = 0.5267, 7 = 1.0263. Its 
linear dependence on x is deceptive, its behaviour near small values of x is actually curved 
(short dashed line in Fig. 1(b)) lying below the curve corresponding to P GPOE (x) (see solid 
line in Fig. 1(b)). For < x < 0.5, we have P GPUE (x) ~ 2.5x(l - .95a;). 



V. INTERESTING DICHOTOMIES AND SPECULATIONS ON RH 

Quasi-Hermitian matrix Hamiltonians 

Pseudo-Hermitian Hamiltonians under a definite metric as given below are called quasi- 
Hermitian. Consider 3- and 4-parameter cases of such 2x2 matrix Hamiltonians. 



a (b + ic)/e 
{b — ic)e a 



, H4 = 


a 7/e 




e 


, V = 






_ 7 *e 




1/e 



(18) 



Here a = a + b, (3 = a — 6, 7 = c + id. These matrices, despite being pseudo-Hermit ian 
possess absolutely real eigenvalues. By constructing one and two parameter pseudo-unitary 
transformation matrices and by carrying out the procedure outlined in Sections, 3 and 4, 
the spacing distributions have been obtained. Interestingly, for H 3 the Wigner surmise, 
P (x) has been recovered identically [8] ! Even more interestingly, by defining e = e _K , 
we find a new analytic expression [13] for P(x) that hardly differs from P GUE (x) for k from 
to 0.5. For other values of k the differences between are also not considerable. 

For e = l, the Hamiltonian is Hermitian, by changing e it becomes non-Hermitian such 
that the spacing distribution does not change appreciably. We may therefore take such 
Hamiltonians and interpret them as a smooth perturbation of the Hermitian Hamiltonian. 

We feel that the display of P GUE (x) by a certain class of spectra, though the underlying 
Hamiltonian is not Hermitian (instead it is quasi- Hermitian) is a very remarkable result. 
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FIGURES 




O 1 X 2 3 °-° x °- 5 



FIG. 1. a: Various spacing statistics, P(x), see Eqs. (4,13,17); b: P(x) for < x < 0.5 

Recall that RZs display P GUE (x) so the prospective Hamiltonian is expected to be both 
Hermitian and TRI-breaking. In the light of our dichotomous result, we speculate that the 
Hamiltonian relevant to RH could also be a PT-symmetric (Hamiltonian). Since complex- 
conjugate eigenvalues are found for a PT-symmetric Hamiltonian when PT-symmetry is 
spontaneously broken, so \ ± it n would naturally follow from a PT-symmetric Hamiltonian. 
Furthermore, the vanishing of the norm (PT-norm) of the eigenstates in the domain where 
spontaneous breaking of symmetry occurs can be seen to be directly connected to a crucial 
criterion proposed by Alain Connes [1] for a Hamiltonian which could be relevant to RH. 
Recently, we have found [13] that for RH a few Hermitian Hamiltonians proposed so far [14] 
do not even possess a discrete spectrum. 

VI. CONCLUSION 

The present work tries to answer the basic question as to what could be the random 
matrix theory of currently researched pseudo- Hermitian Hamiltonians. The new "universal- 
ities" in Eqs. (13) and (17) as displayed in Fig. 1 are distinctly different from the usual ones, 
resulting in weaker repulsion among the energy-levels. Like the well established GOE,GUE 
and GSE, here it remains to be proved that the claimed results for 2 x 2 matrices would 
actually stay on at least as good approximants for the N x N case. This is an open challenge 
an answer to this would actually take us from "universalities" to universalities. 

However, we feel that 'the weaker level repulsion for small spacings' is the essence of 
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pseudo-Hermiticity. And since the pseudo-Hermiticity has been re-cast in terms of more 
physical PT-symmetry, an observance of weaker level repulsion at small spacings would 
suggest an individual violation of P and T symmetries and joint invariance of PT. 

PT-symmetric (pseudo-Hermitian) Hamiltonian as a physical model may now be far, 
nevertheless signature such as weaker level repulsion would give way to such Hamiltonians. 

Two interesting dichotomies have also been presented - owing to one of them we have 
speculated that the prospective Hamiltonian for the Riemann Hypothesis would rather be 
PT-symmetric (pseudo-Hermitian). 

Lastly, we would like to point out that in contrast to the Ginibre matrix ensembles and 
some more suggested by Haake [9], we consider such complex matrices (pseudo-Hermitian) 
which give conditionally real eigenvalues and we define level spacing as {E\ — E 2 ) and not 
as modulus of difference of two complex eigenvalues. 
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